Breaking rotational symmetry in two-flavor color superconductors 
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The color superconductivity under flavor asymmetric conditions relevant to the compact star 
phenomenology is studied within the Nambu-Jona-Lasinio model. We focus on the effect of the 
deformation of the Fermi surfaces on the pairing properties and the energy budget of the super- 
conducting state. We find that at finite flavor asymmetries the color superconducting BCS state is 
unstable towards spontaneous quadrupole deformation of the Fermi surfaces of the d and u quarks 
into ellipsoidal form. The ground state of the phase with deformed Fermi surfaces corresponds to a 
superposition of prolate and oblate deformed Fermi ellipsoids of d and u quarks. 



I. INTRODUCTION 

A possible outcome of the de-confinement phase tran- 
sition in hadronic matter, when it is compressed in the 
centers of compact stars to densities several times the nu- 
clear saturation density, is the formation of a two-flavor 
(u and d) quark matter coexisting with leptons under 
(3 stability and global charge neutrality yj. It is be- 
lieved that the strong interaction among quarks in the 
attractive channel(s) pairs the u and d quarks of two 
color (leaving one color unpaired) to form a Bardeen- 
Cooper-Schrieffer (BCS) type color superconductor (the 
so called 2SC phase) 0, LJ In the asymptotic regime 
of high-densities rigorous perturbativc results can be ob- 
tained in the weak coupling regime ; at the intermedi- 
ate densities, relevant for compact stars, one has to rely 
upon effective- field theory models of QCD, such as the 
Nambu-Jona-Lasinio (NJL) model. 

In compact stars the color superconductivity is likely 
to occur at finite isospin asymmetry, where the u and 
d quarks fill two distinct Fermi spheres. The separation 
of their Fermi energies is of the order of the electron 
chemical potential as required by the charge neutrality 
and the stability against weak processes d —> u + e + v 
and u + e — > d + u, where e, v and v refer to elec- 
tron, electron neutrino and anti- neutrino. Large sep- 
arations of the Fermi energies (of the order of pairing 
gap) and corresponding Fermi momenta implies an in- 
coherence of the phase space for the u and d quarks, 
which will eventually destroy the 2SC superconductiv- 
ity. If the occupation of the single particle states of the 
non-interacting system is preserved when the pairing in- 
teraction is switched on, i.e. the Fermi spheres arc filled 
with quasi-particles in an isotropic and homogenous man- 
ner, the phase-space coherence (at finite separations) can 
be maintained by pairing Cooper pairs with finite cen- 
ter of mass momentum [Larkin-Ovchinnikov-Fulde-Ferrel 
(LOFF) phase] 5] . This is in line with the basic prescrip- 
tion of the Landau Fermi liquid theory which states that 
in strongly interacting Fermi liquids the occupation prob- 
abilities of interacting quasi-particles remain the same as 
for the non-interacting system when the interactions are 
switched on adiabatically. 



The loss of the condensation energy due to the sepa- 
ration of the Fermi spheres of the u and d quarks, how- 
ever, might favor a non-Fermi liquid occupation of the 
fermionic states, where the condensation energy is main- 
tained at the cost of extra kinetic energy. One possi- 
bility is the deformation of the Fermi spheres (DFS), 
to the lowest order, in the ellipsoidal form (hereafter - 
DFS phase) 0. While deviations from the spherically 
symmetric form of the Fermi surfaces cost extra kinetic 
energy, the gain in the condensation energy due to an 
increase in the phase-space overlap between the states 
that pair can stabilize the system. Note that we as- 
sume that the pairing interaction is isotropic in space; 
such a deformation spontaneously breaks the rotational 
0(3) symmetry of the original Lagrangian down to 0(2). 
Both LOFF and DFS superconducting phases break the 
global space symmetries. Another possibility that main- 
tains these symmetries is the rearrangement of the Fermi 
surface of the <i-quarks [interior gap pairing (IGP)] |7(. 
The IGP assumes <i-quark holes located in a strip around 
the Fermi surface of the u-quarks are lifted to the excited 
states above the Fermi surface of the d quarks and the 
pairing occurs at the shores of the it-quark Fermi sea • 

A robust feature of the 2SC superconductors is the ap- 
pearance of the crystalline color superconducting (CCS) 
state (the analog of LOFF phase) in a wide range of 
chemical potential asymmetries H U ^| 0, 0, \l3L lldl 
IfH Pill IrH ll"l . In the CCS state the Cooper-pairs carry 
finite center-of-mass momentum, their order parameter 
varies periodically in space and, hence, the CCS sponta- 
neously breaks t he global sp ace symmetries of the original 

bcs state H H EMl ElEl El El 13 El El o ne of 

the consequences of the broken space symmetries is the 
existence of new massless Goldstone modes (phonons) in 
the CCS phase El El El El; 

another consequence is 
that, the glitches in compact stars (rotational anomalies 
seen in the timing data from pulsars) could originate in 
the CCS phase, if the compact stars feature such a phase 

This paper studies numerically the color superconduct- 
ing DFS phase; it extends the previous non-relativistic 
analysis [6| to relativistic systems and interactions spe- 
cific to the color superconductors. The main result of this 
analysis is that the spatially homogenous 2SC supercon- 
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ductor is unstable towards formation of the DFS phase 
with quadrupole deformed Fermi surfaces. Clearly, this 
observation does not identify the true ground state of a 
2SC superconductor. Studies of the interplay between 
different non-BCS phases sketched above are needed to 
answer this question. For applications to compact stars 
these phases need to be studied under (3 equilibrium con- 
ditions (for example see (13. 12(1 l2lj ) . 

The paper is organized as follows. In Sec. 2, start- 
ing from the NJL Lagrangian in the chirally symmetric 
phase, we derive the thermodynamic potential of the fla- 
vor asymmetric 2SC phase at finite temperatures. Self- 
consistent equations for the gap and the partial densities 
of the u and d quarks are obtained. The numerical solu- 
tions of these equation and the resulting phase diagram 
of the color superconducting DFS phase are shown in Sec. 
3. Sec. 4 contains our conclusions. 



II. THERMODYNAMICS OF COLOR 
SUPERCONDUCTING DFS PHASE 

Our treatment is based upon the Nambu-Jona-Lasinio 
(NJL) with two flavors (Nf — 2) and three colors (N c = 
3) [22]. We shall assume from the outset that the color 
superconducting phase emerges in the chirally symmetric 
phase. The Lagrangian density of the model is then given 
by 

(1) 

where C — ij 2 ~f° is the charge conjugation operator, t 2 
is the second component of the Pauli matrix acting in the 



SU(2) / flavor space, Xa is the antisymmetric Gell-Mann 
matrix acting in the SU(3) C color space. The coupling 
constant G\ stands for the four-fermion contact interac- 
tion. For the sake of simplicity we neglect the effects 
of the mass gap in the quark-anti-quark channel below 
the critical line for the chiral phase transition on the 
Cooper pairing. The generalization of the Lagrangian 
(|T|) , which leads to coupled gap equations in the di-quark 
and Cooper channels, is straight forward. The common 
Ansatz for the order parameter in the 2SC phase is 



A cx (^ t (x)Cj 5 t 2 X 2 ^(x)), 



(2) 



where A2 is the second component of the Gell-Mann ma- 
trix. The Ansatz for the order parameter [Eq. J2Jl] im- 
plies that the color SU(3) C symmetry is reduced to SU(2) C 
since only two of the quark colors are involved in the pair- 
ing while the third color remains unpaired. More com- 
plicated Ansaetze would allow for a spin 1 pairing of the 
quarks of the remaining color [23l |24| . 

As is well known, the gap equation and the partial 
densities of the up and down quarks can be found from 
the fixed points of the thermodynamic potential density 
ft: 



dfl on 



(3) 



the flavor index f = u,d refers to up (u) and down (d) 
quarks. In the superconducting state the self-consistent 
solution of equations corresponds to a minimum of 
the thermodynamic potential of the system. For the La- 
grangian density defined by Eq. Q and the pairing chan- 
nel Ansatz Eq. (01, the finite temperature thermodynam- 
ical potential f2 per unit volume is 



d 3 p 1 
(2Vp 2 



Tr In 



q( Su(iuJ n ,P) S^iiuJn,]}) 



A 2 
4Gi' 



(4) 



where [3 is the inverse temperature. The matrix structure of the inverse Matsubara propagators S^^iuJmP), (l,k = 
1, 2) reflects the Nambu-Gor'kov extension of the particle-hole space to account for the pair correlations. The elements 
of the Nambu-Gor'kov matrix are 2x2 matrices defined as diag [S^i] = {$ + Hul°,i> + Md7°), and diag [S22 ] — 
H> - Hul a ,i> ~ Md7°); diag [S^\ = (A75T2A2, A75T2A2) and diag [S^ 1 ] = (-A*7 5 t 2 A 2 , -A*75T2A2) with off-diagonal 
elements zero. Upon carrying out the traces in the spin, flavor and color spaces and the fermionic Matsubara 
summation over the frequencies uj n one finds 



n 



-glog (1 + e 



Eij + -log (1 



A 2 

4Gi' 



(5) 



where the indeces i,j = (+, — ) sum over the four branches of the paired and unpaired quasiparticle spectra defined, 
respectively, as £±± = (p± /j,) ±5/1 and E±± = ^ (p ± n) 2 + | A| 2 ± <5/x, where — (fi u — /id)/2 and \x = (/!„ + Md)/2 
with jj, u and /id being the chemical potentials of the up and down quarks; s +J = 1 and s_j = sgn(p — Expressions 
analogous to JSJ were derived in refs. pH l26t I2H l29|| for finite/zero temperature and flavor symmetric/asymmetric 
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cases. The variations of the thermodynamic potential (0 provide the gap equation 



A = 8G1 
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where f(£%j) are the Fermi distribution functions and the 
upper/lower sign corresponds to the u/e?-quarks. The 
changes in the thermodynamic potential due to the phase 
transition to the superconducting phase at constant (3 
and pj are equivalent to the changes in the free energy 
at constant (3 and pf when the potentials are expressed 
in terms of appropriate thermodynamical variables. The 
free energy F is related to the thermodynamic poten- 
tial SI by the relation F = fl + p u p u + pdPd- Below, 
we choose to minimize the free energy of the system 
at constant temperature and density of the matter as 
a function of the flavor asymmetry parameter defined as 
a = (pd — Pu)/{pd + Pu)- Thus, fixing the net density 
of the quark matter and the flavor asymmetry, fixes the 
chemical potentials p f by means of Eqs. 10 and (J7J. The 
chemical potentials are isotropic if the Fermi surfaces of 
the species are spherical. To accommodate the possibility 
of their deformation we expand the chemical potentials 
in Legendre polynomials with respect to an (arbitrary) 
axis of spontaneous symmetry breaking of the spherical 
symmetry of the Fermi surfaces: 

oo 

M/=E^ P <( Cos0 )- ( g ) 

1=0 

The I = terms are constants which merely renormal- 
ize the chemical potentials from their values for the case 
of vanishing pairing interactions. Since the interactions 
arc translationally invariant, the I — 1 terms can not 
contribute, since they manifestly break the translational 
symmetry. (The breaking of translational symmetry can 
be caused by the kinetic energy contribution, as is the 
case for the LOFF/CCS phases.) We shall keep below 
the leading order I = 2 terms which break the rotational 
symmetry of the system by deforming the Fermi surfaces 
into an ellipsoidal form. The corresponding expansion 
coefficients are treated below as variational parameters 
to minimize the free energy of the color superconductor. 
With the leading order in the deformation terms kept, 
the chemical potentials of the u and d quarks can be cast 



in the simple form 

p f = fi f [l + (e s ±e A )(Co S 9) 2 ] , (9) 

where the lower/upper sign corresponds to up/down 
quarks; e s/A = {3/2)[p 2d / Pod ± M2«/Mo«] and p f = 
Pfo — M/2/2. The parameters es and ea describe the de- 
formations of the Fermi surfaces; the symmetric part (es) 
is the measure of the conformal expansion/contraction of 
the Fermi surfaces; the antisymmetric part (sa) describes 
the relative deformations of the Fermi spheres of u and 
d quarks. Note that the expansion © need not conserve 
the volume of a Fermi sphere; the isotropic expansion 
parameters are self-consistently found form the normal- 
ization to the same matter density as in the case for the 
undeformed state. 



III. THE PHASE DIAGRAM 

As the NJL model is non-renormalizable, the momen- 
tum integrals in the gap equation 10 need to be regular- 
ized by a cut-off; we employ a three-dimensional momen- 
tum space cut-off |p| < A. The phenomenological value 
of the coupling constant G\ in the (qq) Cooper channel 
is related to the coupling constant in the (qq) di-quark 
channel by the relation G\ = N C /(2N C — 2)G; the latter 
coupling constant and the cut-off are fixed by adjust- 
ing the model to the vacuum properties of the system 
[HlH |27|]. We employ the parameter set d = 3.10861 
GeV" 2 and G1A 2 = 1.31 from Ref. Hj. For this param- 
eter set the value of the gap is 40 MeV at the baryonic 
density pb = 0.31 fm~ 3 which corresponds to the chemi- 
cal potential value p = 320 MeV in the flavor symmetric 
state at T = 0. Figure 1 summarizes the main features of 
the DFS phase in the physically relevant regime of flavor 
asymmetries 0.1 < a < 0.3 likely to occur in the charge- 
neutral matter under /3 equilibrium. The color supercon- 
ducting gap (left panel) and the free energy difference be- 
tween the superconducting state and normal state (right 
panel) are shown as a function of deformation parameter 
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sa for several flavor asymmetries at the baryonic den- 
sity p B = 0.31 fm" 3 and temperature T = 2 MeV. The 
conformal deformation has been constrained to £5 = 0. 
The gap is normalized to its value in the flavor symmet- 
ric and undeformed state A o = A(a = 0, ea = 0) = 40 
MeV. The free-energy difference SF — Fs~- Fn is likewise 
normalized to its value in the flavor symmetric and un- 
deformed state SF(a = 0, e A = 0) = 37 MeV fm" 3 . The 
/3 equilibrated quark matter requires an excess of the d 
over u quarks, therefore the range of the flavor asymme- 
try is restricted to the positive values. The deformation 
parameter ea assumes both positive and negative values. 
The main features seen in Fig. 1 are: (i) for a fixed a/0 
and ea > 0, the gap is larger in the DFS state than in 
the ordinary BCS state (ea = 0); (ii) the minimum of the 
free energy corresponds to the DFS state with ea — 0.25, 
and its position weakly depends on the value of a. 




FIG. 1: The color superconducting gap (left panel) and the 
free energy (right panel) as a function of relative deformation 
parameter ea for the values of the flavor asymmetry a = 
(solid lines), 0.1 (dashed lines), 0.2 (short dashed lines) 
and 0.3 (dashed dotted lines) at density pB = 0.31 fm~ 3 and 
temperature T — 2 MeV. 

Positive values of ea imply a prolate deformation of 
the Fermi surface of the d-quarks and an oblate defor- 
mation for the Fermi surface of the u-quark. Figure 2 
illustrates the Fermi surfaces for the parameter choice 
a = 0.3 and ea = 0.25, which corresponds to the de- 
formation at which the minimum in the free energy is 
obtained. 

As can be seen from Fig. 2, the deformation of the 
Fermi surfaces leads to an area in which the Fermi sur- 
faces of u- and d-quarks touch each other, which means 
that for this part of the phase space there is no blocking 
effect reducing the pairing correlation. As a consequence, 
the superconducting gap is larger than for the configura- 
tion with spherical Fermi surfaces leading also to more 
attractive free energy. 




FIG. 2: The deformed Fermi surfaces of the d and u quarks 
for the flavor asymmetry a = 0.3 and deformation parameter 
ea = 0.25. The right part of this figure refers to the Fermi sur- 
face of the d-quarks, exhibiting a prolate deformation, while 
the oblate shape on the left corresponds to the Fermi surface 
of the u- quarks. 

The behavior of the gap as a function of a and 
ea is best understood in terms of the symmetric and 
anti-symmetric combinations of the quasiparticlc spec- 
tra E s = (E-+ + £— )/2 and E A = (£_+ - )/2. 
Since we assume that £5 = the symmetric combina- 
tion of the spectra is weakly affected by the deformation 
[~ 0(eaSp)}. In the case Ea = one recovers the BCS 
result with perfectly overlapping Fermi surfaces. The ef- 
fect of finite Ea at ea = is to induce a phase-space de- 
coherence in the kernel of the gap equation; this blocking 
effect reduces the magnitude of the gap. 

As the chemical potential shift dp, and ea contribute to 
Ea with different signs, switching on finite ea > acts 
to reduce Ea and the magnitude of the gap increases 
due to the restoration of the phase-space overlap. Once 
the condition ea > 5p is satisfied further increase in the 
deformation acts to increase the de-coherence and, thus, 
decrease the magnitude of the gap. 

For negative values of ea the deformation and the shift 
in the chemical potentials contribute to Ea with the same 
sign, therefore the deformation acts to increase the phase- 
space de-coherence which reduces the magnitude of the 
gap. These features are seen in Fig. 3, where we show the 
pairing gap as a function of the flavor asymmetry a and 
the relative deformation of the Fermi surfaces ea for the 
same density and temperature as in Fig. 1. [We assume, 
as for Fig. f , that the conformal expansion/contraction 
of the Fermi surfaces is absent, eg = 0]. 

Figure 4 shows the difference between the free ener- 
gies of the superconducting and normal states 5F nor- 
malized to its value in the flavor-symmetric BCS state 
SF o = 5F(a = 0,£a = 0). Numerically, the con- 
tribution from the entropy difference between the nor- 
mal and superconducting phases is negligible. Small 
ea > 0, a ^ perturbations from the flavor symmet- 
ric/undeformed state increase 5F signaling an instability 
of the BCS state towards deformations of the Fermi sur- 
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FIG. 3: The color superconducting gap as a function of flavor 
asymmetry a and the relative deformation parameter ea at 
density ps = 0.31 fm -3 and temperature T — 2 MeV. The 
pairing gap is in units A(a = 0, ea = 0) = 40 MeV. 




FIG. 4: The difference between the free energies of the super- 
conducting and normal phases as a function of flavor asym- 
metry a and the relative deformation parameter ea- The free 
energy scale is normalized to its value 8F(a = 0, ea — 0) = 37 
MeV fm" 3 . 

faces. 

For flavor asymmetries a > 0.3 and positive values 
of the deformation parameter the superconducting state 
can exists only in the DFS phase. For negative values of 
sa there is no energy gain related to the deformation of 
the Fermi surfaces (except a marginal effect in the limit 
a — > and sa — * 0). Quite generally the shape of the 
free energy surface reflects the functional behavior of the 



gap as a function of the parameters a and ea, which is 
due to the dominant contribution from the condensation 
energy (cf. Fig. 1). 

IV. CONCLUDING REMARKS 

Color superconducting 2SC phase becomes unstable 
towards formation of a superconducting phase with de- 
formed Fermi surfaces of up and down quarks for finite 
separations of their Fermi levels. Small asymmetries in 
the populations of the of the up/down quarks are suffi- 
cient to achieve the bifurcation point where the form of 
the Fermi surfaces changes from the spherical to the el- 
lipsoidal form. Although the departures from the spher- 
ical form of the Fermi surfaces costs extra kinetic en- 
ergy, there is net energy gain because the deformation 
increases the phase space coherence between the fermions 
forming Cooper pairs and, therefore, the gain in the po- 
tential energy. As mentioned in the introduction, the 
present study does not establish the true ground state of 
the 2SC superconductor under flavor asymmetry. This 
can be done by a careful comparison of the condensa- 
tion energies of the DFS and CCS phases. With the 
form of the lattice of the CCS phase established in the 
Ginzburg-Landau regime |13( it becomes feasible to carry 
out a combined study of the DFS and CCS phases [30j . 

For applications to the compact stars, our treatment 
needs an extention which will include the (3 equilibrium 
and charge neutrality among the quarks and leptons. The 
onset of the DFS phase (or its combinations with non- 
BCS superconducting states) will affect the properties of 
the compact stars in a number of ways. One straightfor- 
ward implication is the modification of the specific heat 
of the superconducting phase which affects the thermal 
cooling of compact stars: we anticipate that the spe- 
cific heat will be suppressed linearly in the DFS phase as 
compared to the exponential suppression in the ordinary 
BCS-state. The spontaneous breaking of the rotational 
symmetry implies emergence of massless bosonic modes 
(Goldstone's theorem). The deformation of the Fermi 
surfaces opens new channels for neutrino radiation via 
the bremsstrahlung process of the type u — > u + v + v 
and d — > d + v + v since the phase-space probability for 
a transition of a quasi-particle from one point on the de- 
formed Fermi surface to another (say, from the equator 
to the pole) is non-zero. 
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